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1. SCOPE 


As an integral part of the Large Area Crop Inventory Experiment (LACIE) 
which is sponsored by the Department of Agriculture (USDA) , the National 
Oceanic and Atmospheric Administration (NOAA) and the National Aeronautics 
and Space Administration (NASA) , continued building of daily and monthly 
data bases to test various advanced wheat yield models is an active objec- 
tive and goal. The People’s Republic of China (PRC) is one of the largest 
agricultural countries in the world. Winter and spring wheat is one of its 
major crops; therefore, a knowledge of China’s crop production is very 
important to the LACIE _ Project. The China data base consists of meteorolo- 
gical information in the form of daily precipitation and maximum and mini- 

1 

mum temperatures and their monthly totals and averages, respectively, for 
about 95 weather stations covering the entire winter and spring wheat produc- 
tion area of China over the period 1965 through 1975. The dally and 
monthly meteorological information was built up from synoptic (4 or 8 times 
dally) weather observations available on magnetic tape from the National 
Climatic Center (N.C.C.) in Asheville, North Carolina. 

There are three major sets of computer algorithms for building the 
China data base, that is: 1) the least squares Interpolation; 2) modified 
version of the max-min program; and 3) the AOV procedure. 

The least squares interpolation algorithm is used to generate input 
variables; the modified version of the max-mln program is used to generate 
daily and monthly precipitation and max-mln temperatures; the analysis of 
variance algorithm is used to carry out a simple statistical hypothesis 
test on the digital estimated data. The China data base is considered 
complete and unique. It may be the only one of this kind of meteorolo- 
gical data existing outside, or inside the People’s Republic of China. 
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2 . APPLICABLE DOCUMENTS 


Action Documentation -.(AD) 63-1347-3763-54, China Data Base, 

June 28, 1977 

AD 63-1557-4963-22, YES Data Base Maintenance, Oct 17, 1977. 

Job Order 73-763 

"As-Built” Design Specification for Historical Daily Data Bases 
for Testing Advanced Models (JSC-12891; LEC-10572), April 1977. 
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3. System Description 


CHINA SYNOPTIC DATA SOURCE 

A total of 32 tapes were obtained by the Center for Climatic and 
Environmental Assessment (CCEA) from the National Climatic Center (NCC) 
in Asheville, North Carolina, covering World Meteorological Organization 
(WMO) blocks, 53, 54, 57 and 58 in China for the years 1965 through 1975. 
,Data were contained in 3 tapes per year except for 1968 which has only 
;2 tapes. These tapes, containing synoptic data (i.e,, 4-8 observations 
per day) for all reporting weather stations within these blocks, are 
recorded in NCC deck format 9685. 


CHINA DATA BASE DESIGN 

A total of 95 WMO weather stations providing uniform coverage of the 
entire winter and spring wheat production area of China (see map on 
page 3.3), was selected on the basis of relative reliability and con- 
sistency of record. Daily precipitation and max-mln temperatures are the 
major meteorological variables generated for the 95 representative wea- 
ther stations covering the entire winter and spring wheat production area 
of China over the period from 1965 to 1975. Daily data were stored in the 
magnetic tapes in the following format: 


Character String 
(per record) 

1-6 
7-8 
9-10 
11 - 12 

13 

14 - 16 
17 - 19 
20 - 22 
23 - 26 


Content 


Station number 

Year 

Month 

Day 

Blank 

Daily Maximum temperature* 
Daily Minimum temperature* 
Blank 

Precipitation ** 


* Temperatures are in whole degrees Fahrenheit with a leading minus 

sign if negative; missing temperature observations are denoted by 999. 
** Precipitation is in hundredths of inches with no decimal points; 

missing precipitation is indicated by 9999. The above format is known 
as the Modified Kansas Format (as opposed to NCC Format) . 
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Monthly data were generated by a separate computer program; output is 
available in hard copies containing the following information: 

1. Station 

2. Year 

3. Month 

4. Monthly average maximum temperature 

5. Monthly average minimum temperature 

6. Monthly total precipitation 

7. Total number of observations 

8. Number of missing observations 

9. Total annual precipitation 

This section describes the system software. Appendices A, B, and 
C contain detail flowcharts for subroutines used. 
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3.1 HASDWARE DESCRIPTION 


All computer algorithms described in this document were written in 
standard FORTRAN language and will operate on any IMS 360/370 system 
with FORTRAN compiler and with' 1600 BPI magnetic drive. They will also 
operate on mini-computer systems such as PDF 11/50 with FORTRAN compiler 
and 800 BPI with minor change. 

3.2 SOFTWARE DESCRIPTION 

3.2.1 THE DERIVATION OF INPUT VARIABLES 

The |least sqxaares interpolation computer algorithm was designed to 
generate a complete set of monthly temperature ranges. The monthly 
average daily temperature ranges for each weather station are needed as 
input variables to the max-mln program to generate the daily maximum 
and minimum temperatures for each month. The program consists of one 
main program and seven subroutines i.e., SOM, COR, LSTSQ, MEQSOL, GRAPH, 
STANDV AND MULTR, and one IBM built-in library routine, PLTSCT, which is 
used to plot the scatter diagram for the input data; the above subroutines 
can be found in computer literature and journals. 

SUBROUTINE SUM 

This subroutine is used to calculate the total of input observations such 
as X, X^. and Y and Y^_. 

SUBROUTINE COR 

This subroutine is used to calculate the correlation coefficient (r) of 
the input data, the correlation coefficient is used as a criterion to 
determine whether a straight line or a 2nd degree curve fits a given set 
of data; that is, if jr [> 0.7, a straight line should be fitted; 
otherwise, a higher degree curve is required. 
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SUBROUTINE LSTSQ 


This subroutine is used to determine the best fitting equation by 

establishing the proper matrix equation AX = B 

Tor example, if a linear least squares equation is needed 

then. 



where are the coefficients of the least squares equation. For 
higher order least squares polynomials, the elements of different 
matrices will be changed automatically by this subroutine. 

SUBROUTINE MEQSOL 

This subroutine is used to solve the matrix equation AX = B which has 
been determined by the subroutine LSTSQ. ORIGINAL PAGE IS 

OP POOR QUALiry 

SUBROUTINE GRAPH 

This subroutine is used to plot the least squares equation in the X - Y Plane. 
SUBROUTINE MDITR 

This subroutine is used to calculate the multiple correlation 
coefficient (RMULT) which serves as the criterion to determine if a 2nd 
degree or higher order least squares equation is needed. 

SUBROUTINE STANDV 

This subroutine is used to calculate the standard deviation of a given 
set of data points; this subroutine is needed in the subroutine MULTR. 

3.2.2. THE COMPUTATION ALGORITHMS 

There are two major computer algorithms for building the China data 
bases. One for generating the daily data and the other for generating 
monthly data. All programs are written in standard FORTRAN language. 
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The computer program for daily data is a modified version of the max-min 
program which is basic to the non-linear statistical model developed by a 
research team at Kansas State University led by Dr. Arlin Feyerherm, 
Professor of Statistics. Complete information about the max-min program 
can be found in the "As-Built Design Specification for Historical Daily 
Data Bases for Testing Advanced Models'*, prepared by Buddy H. Jeun and K. 
Williams and published by the National Aeronautics and Space Administration 
(NASA), Lyndon B. Johnson Space Center, Houston, Texas 77058, April 1977. 

i 

In generating dally precipitation for the PRC, the precipitation amounts 
reported at synoptic hours 0600, 0900, 1200, and 1800 (GMT) were totaled. 
(When a similar application was carried out for the USSR only 0600 and 1800 
(GMT) observations were summed) . The reporting conventions differ region- 
to-region and must be verified locally by summing the data by day, month, 
and year. The computer algorithm that generates the monthly meteorological 
data, i.e. , average maximum and minimum temperature and total precipitationi 
consists of one main routine and five subroutines, which are SORTMD, DMEAN, 
TOTAL, ANP and PBAR, and one built-in function, i.e., EXTERNAL SIGN* 
which is used to distinguish if the observation is blank or zero, since 
blank denotes a missing observation for maximum and minimum temperatures 
and for precipitation, the missing obseirvation is denoted by 9999. 

' SUBROUTINE SORTMD 

This subroutine is used to sort the daily observations by month; in 
doing so, for every station, the maximum and minimum of temperatures and 
precipitation for every month are sorted out. 

SUBROUTINE DMEAN 

This subroutine is used to calculate the average maximum and minimum 
temperatures for each month. At the same time, it also keeps track of the 
number of missing temperatures by the argument ICOUNT and KOUNT. In this 
subroutine, ICOUNT represents the number of missing observations for the 
maximum temperature and KOUNT represents the number of missing observations 
for the minimum temperature. 
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SUBROUTINE TOTAL 


This subroutine is used to calculate the sum of a given set of observations. 
Input arguments are N, D2 and output arguments are SUM and KNT, where N is the 
number of observations. D2 is an array of N observations. SUM is the total 
value of N observations, KNT is the number of missing observations. 

SUBROUTINE PBAR 

This subroutine is used to calculate the monthly total of precipitation to 
each station. Input argument PC is an array of N observations of pre- 
cipitation. Output argument PAV is the monthly total of precipitation. 

SUBROUTINE ANP 

This subroutine is used to calculate the annual precipitation to each weather 
station. Input argttments are STATN, YEAR, N and VP; output arguments are BP 
and the missing observations are represented by ICOUNT. 

3.2.3. THE EVALUATION OF OUTPUT VALIDITY ^ 

OF POOR QUALUi 

The analysis of variance. 

This analysis of variance computer algorithm was designed to test the 
reliability of the monthly data derived from the modified version of the 
max-min program. A few weather stations were selected for a study comparing 
the estimated means with PRC published means obtained from the USAF through 
the National Climatic Center at Asheville, North Carolina and other reliable 
records from the Hong Kong Royal Observatory and Chinese educational institu- 
tions obtained through personal correspondence. The six stations for which 
comparative monthly temperature data are available are shown on the map on 
page 3-9. They are confined to the more humid eastern portions of China but 
they do give a good temperature range north-south across the major wheat- 
producing area. The tables and graphs (page3-10 through 3-21)include individual 
monthly average temperatures (EST) as summarized from the computed daily 
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max-mln values alongside the most reliable check data (PUB) available. 
Whether the check data are true internal published amounts or "usually 
reliable" figures is not known. An analysis of variance (page 3-22 
through 3-24) indicates that the estimated mean monthly temperatures are 
not significantly different from the available check data. The null hypo- 
thesis of the analysis of variance is tenable. 

For precipitation, only four stations were supplied for comparison (pages 
3—25 through 3—29) . The WMO station numbers are the same as in dally 
source tapes but it Is suspected the pairs of records compared may come 
from different, but proximate, rain gauges. Nevertheless, the monthly 
and annual precipitation totals compare favorably. 
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SELEtTED WMO WEATHER STATIONS, 






MONTHLY AVERAGE TEMPERATURE 


CHI-NAN 

#54823 



1965 

1966 

1967 

1968 

1969 

1970 

MONTH 

EST 

PUB 

EST 

PUB 

EST 

PUB 

EST 

PUB 

EST 

PUB 

EST 

PUB j 

- \ 

JAN 

-0.38 

-0.1 

-0.63 

0.3 

-2.70 

Mj 

-3.27 

■ 

-4.66 


-3.86 

1 

-1.8 

FEB 

2.70 

2.7 

3.56 

4.2 

-0.18 

0.4 

-3.16 

H 

-4.86 


0.35 

2.7 

MAR 

8.20 

7.8 

6.54 

7.6 

7.25 

8.6 

7.43 

9.4 

3.37 


3.21 

5.2 

APR 

11.88 

12.7 

13.84 

15.1 

14.50 

16.1 

14.84 

16.7 

11.81 

14.0 

13.04 

15.1 

MAY 

22.48 

22.9 

21.88 

22.7 

23.58 

25.5 

21.29 

22.3 

19.68 

21,4 

19.42 

21.5 

JUN 

26.00 

27.1 

24.50 

26.4 

25.60 

27.3 

26.83 

28.4 

23.76 

26.0 

22.80 

24.8 

JUL 

27.23 

27.7 

27.15 

27.8 

27.48 

27.8 

25.56 

28.7 

25.60 

27.5 

24.82 

26.4 

AUG 

24.89 

25.5 

27.15 

28.1 

27.48 

28.1 

25.56 

26.5 

25.60 

26.9 

24.82 

25.9' 

SEP 

23.50 

23.3 

19.47 

20.6 

18.98 

20.2 

21.46 

22.8 

19.92 

21.2 

19.42 

21.2 

OCT 

20.12 

17.6 

14.79 

15.9 

14.66 

16.2 

12.69 

14.2 

14.58 

16.4 

14.31 

16.3 

NOV 

9.14 

9.0 

8.33 

9.4 

4.46 

6.1 

7.27 

8.8 

4.84 

6.7 

5.50 

7.2 

DEC 

i 

1.20 

1.2 

-2.10 

-1.0 

-4.16 

-3.0 

1.42 

3.2 

-1.66 

0.3 

0.46 

1.5 
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MONTHLY AVERAGE 


TA-LIEN 
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1965 

1966 

1967 

MONTH 

EST 

PUB 

EST 

PUB 

EST 

PUB 

JAN 

-5.60 


-6.20 

-4.5 

-6.80 

m 

FEB 

-2.96 


-3.58 

-1.6 

-4.95 


MAR 

1.76 


-0.33 

1.2 

1.52 

3.4 

APR 

7.03 


5.87 

7.9 

7.00 

9-0 

MAY 

14.83 

15.7 

14.26 

15.9 

14.29 

16.6 

JUN 

19.76 

20.1 

17.48 

18.0 

17.40 

19.4 

JUL 

22.00 

23.1 

21.71 

22.6 

21.90 

23.5 

AUG 

22.13 

22.6 

23.34 

24.6 

24.00 

25.3 

SEP 

20.69 

20.7 

17.27 

18.5 

17.11 

18.7 

OCT 

16.19 

15.5 

13.23 

14.2 

13.36 

14.7 

NOV 

5.53 

6.0 

5.47 

6.9 

2.38 

3.9 

DEC 

-2.77 

-2.1 

-5.50 

-4.0 

-6.80 

-5.4 


> 3 -] 




-5.i6 

-5.46 

4.29 

10.81 


15.10 


20.34 

22.26 

22.60 


1969 

EST 

PUB 

-8.90 

m 

-9.28 

D 

-1.82 

-0.2 

6.35 

8.4 

13.02 

15.0 

17.44 

19.2 

21.37 

22.5 

22.34 

24.0 

18.75 

19.8 

S. 


11.83 

14.0 



, 5.72 ■ 7.1 

1 

1-1.65 0.4 


4.22 

-2.25 
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1965 

1966 

1967 

1968 

1969 

1970 

.‘month 

EST 

PUB 

EST 

PUB 


PUB , 

■ 

PUB 

EST 

PUB 

EST 

PUB 

JAN 

3.41 

3.8 

2.36 

2.9 

0.96 

D 

0.35 

■B 

-0.15 

HI 

-0.68 

0.8 

FEB 

4.37 

5.1 

4.90 

5.9 

2.88 

3.8 

0.76 

1.8 

-1.56 

0.3 

4.02 

5.2 

MAR 

7.58 

7.9 

8.81 

9.7 

8.04 

8.8 

8.73 

9.7 

5.57 

7.1 

4.60 

5.9 

APR 

12.10 

12.7 

13.84 

14.6 

13.33 

14.5 

16,00 

15.1 

13.69 

14.6 

12.89 

14.0 

MAY 

20.31 

20.9 

18.61 

19.8 

20.04 

20.7 

20.78 

19.8 

20.10 

21.3 

18.63 

20.0 

JUN 

23.41 

23.6 

23.51 

24.8 

25.51 

26.0 

25.67 

24.9 

22.10 

23.7 

22.01 

23.3 

JUL 

27.99 

28.4 

27.86 

28.9 

27.36 

28.2 

27.20 

26.8 

25.12 

26.6 

26.55 

27.2 

AUG 

25.26 

25.9 

28.56 

29.6 

29.49 

30.5 

26.43 

27.4 

26.98 

28.4 

27.39 

28.2 

SEP 

21.52 

22.0 


21.1 

21.73 

23.0 

21.98 

23.1 

22.96 

24.1 

21.39 

22.6 

OCT 

16.91 

17.0 

15.76 


14.95 

16.4 

14.55 

15.4 

16.21 

17.1 ' 

16.15 

17.4 

NOV 

11.31 

11.9 

9.18 

10.4 

8.31 

9.8 

11.22 

11.9 

6.79 

8.3 

8.86 

10.0 

DEC 

2.21 

3.4 

2.26 

3.0 

-0.56 

-0.1 



5.92 

B 



3.32 

4.7 
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1965 

1966 

1967 

■ ■ ■■> " ' 1 
1968 

1969 

1970 

MONTH 

EST 

PUB 

EST 

PUB 

EST 

PUB 

EST 

PUB 



EST 

PUB 

-JAN 

-15.49 


-16.87 

-16.9 

-14.27 

-14.4 

-18.43 

-17.8 

-17.72 

-17.8 


-20.5 

1 FEB 

-13.61 



-9.8 

-10.92 

-11.2 

-16.43 

-16.4 

-15.95 

-15.8 


-13.4. 


-4.25 


-5.22 

-4.8 

-1.72 

-1.3 


-1.8 

-4.36 

-4.7 

-8.15 

-8.1* 

APR 

5.22 


6.26 

6.0 

7.48 

7.6 

10.15 

9.5 

6.37 

6.0 

8.10 

! 

7 . 9 ; 

MAY 

1 

15.29 

14.9 

15.78 

16.1 

17.68 

17.4 

15.11 

15.3 

12.93 

13.3 

15.66 

15. 9| 


21.51 

21.3 

20.40 

20.3 

19.00 

18.9 

19.11 

18.9 

18.02 

17.9 

22.39 

21.9 

^9 

22.88 

22.3 

22.52 

22.4 

23.24 

22.9 

24.36 

23.7 

22.75 

22.6 

22.13 

I 

21. 7| 

AUG 

20.83 

20.4 

21.42 

21.5 

23.25 

23.4 

21.01 

21.2 

20.61 

20.5 

22.28 

4 

22 . 1 : 

1 

SEP 

16.25 

15.2 

13.94, 13.8 

15.90 

15.2 

15.07 14.9 

14.44 

14.4 

15.06 

15.2' 

« 

OCT 

8.28 

7.9 

8.38 

8.7 

9.27 

8.7 

6.12' 6.2 

6.65 

7.1 

7.22 

7.5. 

NOV 

-4.30 

-4.1 

-4.08 

-3.2 

-6.73 

} 

-6.1] —3.98 

1 

-3.3 

-6.93 

-6.6 

-2.72 

-2.4 

DEC 

-14.81 

-15.3 

-15.95 

-15.9 

-18.30 

-17.9 

-10.04. -9.2 

-16.81 

-16.1 

-12.42 

-12.6 

1 
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1965 

1966 

1967 

1968 

1969 

1970 

MONTH 

EST 

PUB 

— 

EST 

PUB 

EST 

PUB 

EST 

PUB 

r 

EST 

PUB 

EST 

PUB ■ 
1 


-7.06 

m 

-6.50 

-5.7 

-6.08 

-4.8 

-6.96 

-5.6 

-8.22 

-6.4 

-7,04 

1 

-5.6; 

) 


-1.41 

D 

-2.43 

-2.5 

-4.79 

-3.3 

-6.37 

-4.9 

-7.15 

-5.1 

-4.03 

-2. O’ 


4.50 


1.70 

3.4 

3.63 

4.8 

4.38 

5.7 

-0.15 

2.1 

-0.10 

2.0 

APR 

10.65 

10.9 

10.51 

12.6 

10.45 

12.8 

12.11 

14.2 

10. 3A 

12.9 

11.31 

13.6 

MAY 

20.35 

20.8 

17.79 

19.7 

20.63 

22.1 

17.61 

20.3 

16.96 

19.3 

16.59 

19.7 


23.98 

24.6 

22.28 

23.8 

21.16 

23.7 

23.06 

25.5 

20.92 

23.3 

20.12 

22.8 

JUL 

1 

25.87 

26,9 

24.02 

25.6 

23.63 

25.3 

25.47 

26.9 

.23.13 

24.6 

23.16 

24.8 

1 AUG 

24.19 

24.8 

23.94 

25.0 

23.72 

25.1 

22.68 

CM 

CM 

22.12 

24.2 

22.51 

27.8 

SEP 

20.83 

20.2 

17.35 

- 

18.7 

16.67 

18.3 

17.99 

19.7 

17.56 

19.6 

17.69 

18.9 

OCT 

16.20. 

14.1 

11.42 

12,8 

11.44 

13.0 

9.83 

11.2 

10.19 

12.3 

11.31 

13.0 

' NOV 

5.42 

4.6 

2.50 

3.7 

0.49 

2,0 

2,86 

3.9 

1.38 

2.8 

1.39 

3.0 

1 

1 DEC 

i 

-4.10 

-3.9 

-5.70 

B 

-7.59 

-6.0 

-3.25 

-1.7 

-5.09 

-3.4 

-4.40 

-3.1 


3-18 
























0c-'-!ai'!h"W(''i Mj I tu'riHiti Ni « Yii. I Sd(f|IE 19 X t9 F9 THE CtKIIIUKt UIJI4-0I 



3-19 









MONTHLY AVERAGE TEMPERATURE 


SHANGHAI 

#58367 


t 

1965 

1966 

1967 

1968 

1969 

1970 

MONTH 

EST 

PUB 

EST 

PUB 

EST 

PUB 

EST 

PUB 

EST 

PUB 

EST 

PUB 

1 

1 

JAN 


B 

3.62 

4.3 

1.42 

2.0 

1.85 

. 2.8 

2.10 

3.6 

0.81 

■ ■ 1 

2.0 

’ FEB 



5.55 

6.6 

3.50 

4.0 

0.3^ 

1.8- 

0.93 

2.6 

4.80 

6.0 1 



B 

9.25 

10.2 

8.08 

8.7 

7.94 

8.8 

5.37 

6.8 

4.87 

6.0 1 

APR 


12.2 

13.36 

13.8 

12.88 

13.6 

14.11 

13.7 

12.91 

14.2 

11.67 

13.1 ‘ 
1 

MAY 

19.06 

19.5 

17.64 

19.0 

18.18 

19.3 

18.98 

18.7 

18.55 

19.6 

17.70 

18.7 • 


21.29 

22.2 

22.49 

23.6 

23.91 

24.3 

22.90 

22.6 

21.31 

22.5 

21.19 

22.1 


27.28 

28.1 

26.09 

27.1 

27.03 

28.0 

27.07 

26.7 

25.68 

26.7 

26.73 

27.1 , 

AUG 

26.46 

26.9 

27.44 

28.6 

28.53 

29.5 

26.37 

2 7.. 5 

27.10 

28.0 

26.89 

28.1 

SEP 


22.7 

21.10 

22.3 

23.14 

24,0 

22.57 

23.7 

23.96 

25.0 

23.85 

24.7 

OCT 

18.09 

18.0 

16.66 

18.1 

17.10 

18.4 

16.08 

17.4 

16.42 

18.0 

17.55 

18.5 

a 

NOV 

13.31 

13.9 

11.13 

12.6 

11.04 

12.5 

12.19 

13.3 

9.26 

10.3 

11.22 

11.8 

DEC 

4.26 

4.9 

5.19 

6.1 

-0.04 

1.0 

7.88 

9.5 

3.93 

4.9 

5.87 

6.9 
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ANALYSIS- OF VARIAHCE 


SOURCE 

OF 

VARIATIONS 


£UKI 

^Bams 

MEAH SPHE 

F -RATIO 

BErWEEM GROUPS 

n^mi 

JBSS 

BMS « BS^p-i) 


mm ftRouPS 

N-P 

WSS 

WMS=WSS/^-p) 

p=m/ 
r /WMS 

TOTAL 

N -1 

TSS 




IMSM P 15 THE M OF OTpS 

M IS the « OF OBSERVATifiHS 

•Ri is WE « of fi:fcsaii'ATHN per gj^oup 


TO 

% "* 

WS5 -^(?x,iA 
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ORIGINAL PAGE IS 
OE POOR QUALITY 








ORIGINAL PAGE IS 
OP POOR QUALIIY 


U) 

I 

N3 

U3 


LOCKHEED 


SUM = 


1815.79 


TA-LlEN{#b^662) 

U.S.AIR FORCE HONG KONG ROYAL 


21.80 

23.34 

25.000 

26.67 

• 27.86 

25.000 

35.18 

36.14 

36.000 

45 .60 

45 .86 

48.000 

58.69 

60.26 

59.000 

67.58 

68.18 

68.000 

71.61 

73.58 

73.000 

71.48 

72.68 

76.000 

69,24 

69.24 

69.000 

61.15 

59.90 

58,000 

41.95 

42.80 

43.000 

27.00 

28.00 

30.000 

3.06 




SUM OF GROUP ( 1) IS 597.95 
SUM OF GROUP ( 2> IS 607.84 
SUM OF GROUP( 3) IS 610.00 


SOURCE OF VARIATION 


ANALYSIS or VARIANCE 
DEGREE OF SUM OF 

FREEDOM SQUARE 


MEAN SQUARE 


HETWEEN GROUP 
WITHIN GROUPS 
TOTAL 



6.938 3.A69 
11820.375 358.193 
11827.313 


F RATIO 

0.010 


SINCE FCAL<FTAB»S0 ACCEPT THE HYPOTHESE 
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LOCKHEED 


19.28 

29.44 

40.45 

51.17 
6a. 64 

75.17 
78.58 
75,54 
70,50 

61.18 
41.77 
24.65 


SUM= 


1909.19 

SUMSQ= 

116026.31 

SUM 

OF 

GROUP! 

1) 

IS 

636 . 37 

SUM 

OF 

GROUP ( 

2) 

IS 

636.82 

SUM 

OF 

GROUP! 

3) 

IS 

636 .00 


PEKING(»54511) 

U.S.AIR FORCE HONG KONG ROVAL 


23.18 

23.000 

28.94 

29,000 

39.92 

41.000 

51,62 

57,000 

69 • 44 

68,000 

76,28 

76.000 

80,42 

79,000 

76.64 

76.000 

68 . 3b 

68.000 

57.36 

54.000 

40,28 

38.000 

24.38 

27,000 


SOURCE OF VARIATION 


analysis of variance 

DEGREE OF SUM OF 

FREEDOM SQUARE MEAN SQUARE 


F RATIO 


BETWEEN GROUP 2 0.125 0.063 

WITHIN GROUPS 33 14776,313 447.767 0.000 

TOTAL 35 14776,438 


SINCE FCAL<FTAB,S0 ACCEPT THE HYPOTHESE 


P o 

II 



MONTHLY PRECIPITATION 


(In Millimeters) 



1967 

1968 

1969 

1970 


PEKING 

CHI-NAN 

TA-LIEN 

NANKING 


#54511 

#54823 

#54662 

#58238 

MONTH 



EST 

PUB 

EST 

PUB 

EST 

PUB 

EST 

PUB 

JAN 

2.0 

2.8 

5.8 

7.3 

0.0 

0.7 


9.3 

FEB 

1.7 

2.4 

1.7 

1.5 

5.3 

3.4 


84.1 


10.6 

7.7 

5.8 

3.4 

27.1 

16.0 

50.8 

66.6 

APR 

29.7 

26.8 

41.4 

45.1 

80.2 

70.3 

115.8 

120.0 

MAY 

30.2 

17.4 , 

19.3 

22.1 

42.4 

32.0 

132.0 

155.5 

JUN 

104.6 

129.7 

10.4 

10.0 

26.1 

51.4 

87.1 

91.5 

JUL 

149.1 

173.5 

34.5 

27.9 

105.4 

87.8 

176.2 

218.2 

AUG 

177.1 

180.4 

25.1 

61.3 

101.1 

137.0 

198.8 

155.2 

SEP 

27.4 

26.0 

8.8 

10.2 

15.7 

38.7 

119.3 

154.7 

OCT 

37.8 

19.7 

64.5' 

68.9 

10.4 

7.9 

74.4 

71.1 

NOV 

7.5 

7.0 

37.5 

37.5 

15.4 

32.1 

12.9 

13.7 

DEC 

0.0 

0.0 

40.6 

25.5 

5.6 

4.6 

26.4 

14.7 


577.7 

593.4 

295.4 

320.7 

434.7 

481.9 

1075,4 

1154. 6 1 
j 

TOTAL (In) 

22.7 

23.3 

11.6 

12.6 

17.1 1 18.0 

42.0 

1 

45.0 , 
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4 . OPERATIONS 


4.1 LEAST SQUARES INTERPOLATION 

This program was designed to work on more than one data set which is 
denoted by a parameter NS. 


Each data set contains XP and YP up to 100 observations as xnput data ; 
in each card, punch XP in columns 2-6 and YP in columns 7-11, with one 
decimal point at column 5 and column 10 for XP and YP; respectively. 


At the end of each data set punch a trailer card with -99.0000.0 in 
columns 2-11. 


Example of a program setup; 

// JOB CARD 

//SS EXEC FCALCMGO, PARM.FORM = ’MAP* 
//FORT*SYSIN DD * 

SOURCE PROGRAM OF LEAST SQUARES INTERPOLATION 
//GO-SYSIN DD * 

data cards 

/* 

// 


4.2 THE MODIFIED MAX-MIN PROGRAl-1 

The dally estimates of maximum and minimum air temperatures and the sum 
of precipitation measurements reported each day are obtained according 
to procedures described in "As-Built Design Specifications for Historical 
Daily Data Bases for Testing Advanced Itodels", AD 63-1347-4963-08 
(JSC-12891, LEC-10572, April 1977). 
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4.3 THE MONTHLY DATA BASES ALGORITHM 

This program needs no Input data from cards; all Its Input variables 
are read in from the magnetic tape. 

Daily data are stored in tapes by years, each file contains one year's 
daily data, each of 85-95 stations. 

Example of a program setup: 

// JOB CARD 

/* SETUP TAPE = 1 

//SS EXEC FTG1CL6, PARM.FORT = 'MAP' 

// FORT-SYSIN DD * 

SOURCE PROGRAM OF THE MONTHLY DATA BASES ALGORITHM 
//LKED-SYSIN DD . * 

INaUDE BMDSUB (BMD2992) 

//LKED»BMDSUB .DD DSN^SYS2*BMD*LOAD, DISP = SHR 

//GO»FT 09F001 DD DSN=CH68A, 

// UNIT=TAPE, VOL =(, RETAIN, SER = X06121) , 

// DCB =(RECFM=FB, LBECL=26, BLKSIZE=7800, DEN=4) , 

// LABEL =(44,,, IN), 

// DISP=OLD, SHR) 

//GO-FTlOFOOl DD DSN=CH68B 

// UNIT = TAPE, VOL= SER = X06121, 

// DCB= (RECFM = FB, LBECL = 26, BLKSIgE=7800, DEN=4), 

// LABEL = (45,,, IN) 

// DISP = (OLD, KEKP) 

//GO^SYSIN DD 

/* 

// 
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4.4 THE ANALYSIS OF VARIANCE ALGORITHM 


Twelve input cards (one for each month) are required for running the 
program. Each card contains three estimates of the monthly mean tempera- 
ture based on the KSU algorithm output, the "published" values obtained 
through NCC from ETAC for years available and period of record means 
obtained from various Chinese publications. These values are located as 
follows : 

Colimins 1-5 contains KSU estimates 
" 7-11 " ETAC record 

" 13-18 " period of record means 

and each field contains two digits, i.e., one decimal point punched in 
columns 3, 9, and 15, for KSU, ETAC and period of record means. 

Example of a program setup: 

// JOB CARD 

//SS EXEC FTGICLG, PAEM . FORT — 'MAP' 

//FORT*SYSIN DD’ * 

SOURCE PROGRAM OF AOV ALGORITHM 
//GOoSYSIN DD * 

DATA CARD 

/* 

// 
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FOPT'^AN IV G1 RELEASF 2.0 " MAIN 


DATE = 77321 00/54/57 


PAG! 




0001 

0002 

0003 

0004 

0005 

0006 

0007 

C008 
0009 
0010 - 
001 1 
00 12 
0013 
001 4 

0015 

0016 
00 17 
0018 

0019 

0020 
0021 
0022 

0023 

0024 


C * LEAST SO INTEPPOLATIOH' * 

C w BY 

C * BUDDY H JEUN * 

C j. 4t* « 

c 

'C IN THIS PROGRAM, ALL YOU HAVE TO DO IS TO FIT A SET OF DATA 

C- PTS TO THE COMPUTER.AND THE COMPUTER WILL RETURN A BEST FIT -.EAST 

C SQ EQUATION FOR YOUR DATA PTS, AND THAT RESULT I NG EQUATION IS 

C INDEED THE BEST FIT IN THE STATISTICAL SENSE. 

C AND THE RESULTING EQUATION IS USE FOP INTERPOLATION' 

C 

C THIS PROGRAM CONSISTING OF ONE MAIN PROGRAM AND SEVEN SUBROUTINE 

C SURPROGPAMS.THAT IS SUBROUTINE SUM ,COR, LSTSQ, MEQSOL ♦ AND GRAPH 

C AND STANDV. MULTR ' 

C SUBROUTINE SUM IS USE TO FIND ALL SUMS FOR X AND Y 

C SUBROUTINE COR IS USE TO FIND THE COR. COEFF 

C SUBROUTINE LSTSQ IS USE TO FIND THE N DEG EQUATION 

C SUBROUTINE MEQSOL IS USE TO SOLVE THE N4N SYSTEM OF EQUATION 

C SUBROUTINE GRAPH IS USE TO OUTPUT THE DATA IN THE FORM OF A 

C STRAIGHT LINE OR A CURVE 

C SUBROUTINE STANDV IS USE TO CAUCULAJTE THE STANDARD DEVIATION OF 

C A GIVEN SET OF' DATA PTS 

C SUBROUTINE RMULT IS USE TO FIND THE MULTIPLE CORRALATION COEFF 

C THE CRITERIA TO DETERMINE WHICH ONE IS THE BEST FIT IS ACCORDING 

C TO THE VALUE OF THE MULTIPLE CORRELATION COEFF 

C THE HIGHEST VALUE OF MULTIPLE R IS CONSIDER AS THE BEST FIT 

C 
C 

. . .REAL XP(100),yP( 1001, XL( 2001 ,YU{ 200) ,R=S< 100) 

COMMON/AREAl/SX,SY,SXSQ, SYSQ , SXY/ ARE A 2/ C { 25 ) . A ( 25, 25),B<25) 

REAL XSQ( 1 00) • XCU8( 1 00) 

NS=20 

C NS IS THE NUMBER OF THE DATA SET TO BE FIT INTO THE REGRESSION LINE 

K=1 

9 WRITE<6,20)K 
C HEADING FOR INPUT DATA 

20 FORMAT! • 1 • ,52X ,* THE INPUT DAT #*.r3, • I S 1 * , /60X , • XP • . 1 OX. * YP * // 
,59X,3( •-• ) .8X.3( • -• ) ) 

DO 111 1=1 . 100 
PEAD(5»3) XPU ) ,YP< I ) 

3 FORMAT! 1X.2F5.1 ) 

IF (XP( I ) .FQ.-99 .0 ) GO TO 222 
WRITE! 6,66) XP( I ) , YP( I ) 

66 F0RMAT!56X,F7.2,5X,F7.2) 

111 CONTINUE 
222 NE=I-1 

DO 50 1=1, NE 
XSQ! I ) =XP( I ) 

. XCUB! I )=XP( I )*-k3 

50 CONTINUE 

EN=FLOAT!NE) 

CALL SUM!NE,XP,YP) 

WRI TE!6,7) 

7 FORMAT(///20X, ‘THE RESULT OF CALCULATIONS ARE J • /• +« ,20X .27 ! • •)/) 

WRITE(6,6) NE 
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FORTRAN TV 61 RELEASE 2,0 


MAIN 


DATE = 77321 00/54/57 PAGt 


0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 
004 1 
0042 
004 3 

0044 

0045 

0046 

0047 

0048 

0049 

0050 


0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 


6 F0RMAT(20X, »N=» .IS) 

WRITER 6 . 8 ) SX.SY.SXSO. SVSO.SXY 

8 FORMAT! 20 X. 'BXs" .Fl0i2/20X. *SY= *.F10. 2/2CX, * SXSQ=* .F10.2/20X. 
SYSQ=» .Fl0.2/20X.'SXY=‘.F10o2) 

CALL COP{NE»R) 

WR1TE(6.77 )R 

77 FORMAT! //BOX. ‘THE COR. COEFF IS R=*«F5.2//) 

IF(AaS(R ) .LT .0.7 ) GO TO 444 

C SINCE X AND Y ARE HIGHLY CORRELATED. SO FIT A STRAIGHT LINE 

ANEW=SXSO*SY-SX*SXr 
BNEW=EN*SXY-SX*SY 
DENE W=E N*SXSQ~SX'j'*2 
A1=ANEW/DENEW 
81=BNEW/0ENEW 

WRITE{6.33|A1 .81 - - - 

33 FORMAT(//50X. *THE LINEAR REG LINE I S ? V50X. • Y= ‘ » F8 . 2t 2X . » 4 • 8. 2 . 

+»*X*//) 

WR1TE(6.21) 

‘ XM1N=1.0 

XMAX=26 
DALTAX=1 .0 
X^=XMI N 

666 Y=Al+ai4X 

WRITE (6.22 ) X »Y 
CALL GRAPH! Y) 

X=X+DALTAX 

IF! X-XMAX ) 6 66 .666.317 

C FIND A HIGHF.R DEGREE EQUATION BY THE LEAST SO METHOD 

444 NC=3 

CA LL PLTSC T ( XP . YP . XL. YL . NE » I 00 * 2 . 

+ 36HLEAST SO INTERPOLATION PLOT 
#,36HTHE months 

' +.36HTHE MONTHLY TEMP RANGE 

7 ) 

CALL LSTSQ! XP .YP .NE »NC> 

CALL MULTR! NE .XP. YP.RXV) 

CALL MULTRINE.XSQ.YP. RXSQY ) 

. Cl=C(l) 

C2=C(2) 

C3=C< 3) 

C AFER CALL THE 2ND DEGREE POLYN 

CALL STANDVINE.XP.SDX ) 

CALL STANDV! NE . YP. SOY) 

CALL STANDVCNE.XSQ.SOXSQ) 

C B!I)=C<I)*SD(I )/SDY 

B1=!C2*SDX)/SDY 
B2=!C3+SDXSa)/SDY 

C THE MULTIPLE R FOR 2ND DEGREE IS: 

S2=AaS(RXY*81+RXSQY *18 2) 

R2=S0RT(S2‘) 

NC«4 

CALL LSTSQIXP.YP.NE.NC) 

CALL MULTR! NE .XCUB.YP.RXCUBY) 

C AFTER CALL THE 3RD DE GR E POLYN? 

CALL STANDVINE.XCUB.SOXCUB) 

Bl5s(C!2)*SDX)/SDY 
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FORTRAN IV Gl RELEASE 2»0 


MAIN 


DATE = 77321 


00/54/57 


0069 

0070 

0071 

0072 

0073 

0074 


0075 

0076 


0077 

0078 

25 

0079 

1 1 

0080 

> ■ 

0081 

0062 

„ _ 12 
C 

0083 

0 094 
0005 
0066 
0087 

21 

0088 

201 

0089 

• 

0090 


0091 

0092 

0093 

22 

0094 

0095 

202 

0096 

500 

0097 

319 

0098 

0099 

0 too 
0 101 


0102 

0103 

404 

0104 
0 1 05 

78 

0106 


0 107 

505 

0108 

r 

0109 

0110 
0 n 1 

55 

0 112 

317 

0 113 

0114 

0115 

318 


B2* ( C ( 3 ) *SD XSa ) /SDV 

B3s=<CI4 )»«SDXCU8>/S0y 

THE MULTIPLE P FOR THE 30EGREE IS! 

S3«(RXV4ai+RXSQY*B2+RXCUBV*03 ) 

R3=SQRT(S3) 

IF(R2.GT.R3) GO TO 500 

CALL PLTSCT(XPiYP,XL»YL»Ne, 10 0 f 3, 

+ 36HLEAST SO INTERPOLATION PLOT 
#.36HTHE MONTHS 
+.36HTHE MONTHLY TEMP RANGE 
+ ) 

CALL LSTSQ( XP. YP.NE.NC) 

WRITE (6,25) (I,C(I),1=1,NC) 

FOPMAT(5X, *C( • . I 1 , • ) =»,F15.5) 

WRITE(6 ,11) 

FORMAT! //SOX, *THE LEAST SO EQUATION IS 1 • // ) 
WRITE(6, 12)(C( I). I=1,NC) 


FORMAT(50X, * Ys* 
HEAD ING 
WRITE (6,21) 


FORMAT! " 1 • ,45X •* THE GRAPH 
X» ,5X, 'Y* ,2QX,7a( •_• ) 


F7.2 ,* + * ,F7 .2 , 'X« ' ,F7.2, ♦ + »,F7,2, •X*43* ) 

FIT least SO EQUATION IS!*/ 


» ,2X,6(*_* ) ) 


OF THE BEST 
, 'Y-AXIS'/ 


+4X , 

+ ' +' ,6 ( ‘ ♦ 

XM IN= I .0 
XMAX=26 
OALTA X=1 ,0 
X=XMIN 

Y=C(1 )+C(2)>rX + C(3)<'X**2+C(4)’*X>«‘*3 
CALL GRAPH! Y) 

WRITECf ,22 )X,Y 
FORMAT! IX,F6.2,F8.2) 

X=X+DALTAX 

IF(X-XMAX)20l,20l ,202 
WRITE! 6. 310) 

GO TO 505 

WRITE! 6,319)CI ,C2,C3 

FORMAT (SOX, • Y= • , F7 .2, ' ♦ ♦ , F7 , 2 • ' X* , ' + ' , F 7 . 2 , • X* ♦2 • ) 

XMIN=0.0 

XMAX=26.0 

DALTAXsl .0 

XsXMIN 

Y=C!+C2*X+C3*X+*2 
WRITE(6*78)X,Y 
FORMAT! 60X»2F1 0.2) 

X=X+DALTAX 

IF! X-XMAX)4 04 ,404 ,505 
WRITE!6,55)R2,R3 
FOR'MAT! //2CX. 'THE 
W20X,'THE MULTIPLE 
IF!K.E0.NS) GO TO 
K=K+l 
GO TO 9 
WRITE(6 ,31 e 
F0RMAT(38X, 

IF! K .EO.NS ) 

K=X^1 


Q Q 
tr5 ^ 

RS 

I® 


MULTIPLE P FOR THE 2ND DEGREE ^OLVN I S', F 10. 2/ 
R FOR THE 3RD DEGREE POLVN IS* ,F1 0 .2// ) 

690 


) 

• X-AXIS* ) 
GO TO 690 


PAG 
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GO TO 9 
690 STOP 
END 
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V G1 ~ RELEASE 2*0 SUM DATE = 77321” 00/54/37 

SUBROUTINE SUM(NEiX.Y) 

COMMON/AREAl/SX »SY tSXSQf SYSQ» SXY 
C 

C IN THIS SUBROUTINE SUM, THE INPUT DATA ARE NEtX.Y AND THE OUTPUT 

C ARGUMENTS ARE SX * SY *S XS Q»S YSQ. SXY ; A-L VALUES FOR THE SUMS ARE 

... _C calculated in THIS SUBROUTINE. . . ... 

C 

DIMENSION X(1 00) ,Y(100 1 

SX=0.0 

SY=0,0 

SXSQ=0. 0 

SV£Q = o .0 

SXY=0.0 

DO 10 Isl.NE 

SX=SX + X ( I ) 

SY=SY+Y(I ) 

SXSQ=SXSQ + X( I )**2 
SYSQ=SYSQ+Y ( I 
SXY=SXV+X{ I ) +Y( I ) 

10 CONTINUE 
RETURN 
end 
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000 1 
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0003 

0004 

0005 

0006 

0007 

0008 
0009 
001 0 
001 1 
0012 
00 1 3 

0014 

0015 

0016 
^ 0017 

M 

LO 


SUBROUTINE LSTSQ ( X , Y . N » M) 

C0MM0N/AREA2/C(25) , A( 25. 25),8( 25> 

C 

C LEAST SQUARE CURVE FITTING 

C TO FIT N GIVEN POINTS X(N),Y<N) 8V A POLYNOMIAL 0= 

C M-1 DEGREE, SUM C ( I »7X**t I -1 ) FOR 1=1 TO I=M 

C C{M) IS THE OUTPUT COEFF ARRAY 

C NEEDS SUBROUTINE MEQSOL TO SOLVE FOR AC=B 

C 

REAL X< 100), Y( 100) 

DO 7 1=1 . M 
B{ I )=0, 0 
00 5 K=1,N 

IFUX(K).EQ.O.O).AND.(I.EQ,l))X(K) = 1.0 
5 BU ) = B( 1 )+XCK )♦♦( I-l)*Y(K > 

DO 7 J=1 , M • 

A< I , J )=0, 0 
DO 7 K=l,N 

IFUX(K)^EO.0.0).AND.{(I-fJ).EQ.2))X(K)=1.0 
All , J )=AI I , J )+X(K ) *★( I+J-2) 

7 CONTINUE 

CALL MEQSOL(M) 

RETURN 

END 


SI 
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M^QSOL date a 77321" " ' 00/34/57 


*r 

to 

00 


000 1 

Q002 

c 

c 

c 

c' 

c 

c 

0003 

0004 

0005 

0006 

0007 

0008 


0009 I 

0010 

0011 2 

0012 

0013 3 

0014 4 

0015 

0016 

0017 

0018 • 5 

0019 - 

0020 
0021 

0022 ' 6 

0023 

0024 7 

0025 8 

0026 

0027 _ . 9 

0028 

0029 95 

0030 

0031 ’ 98 

0032 99 

0033 

0034 


SUBROUTINE MEQSOLIN) 

C0MM0N/AREA2/X { 25 ) t A( 25 . 25 ) t B ( 25 ) 

THIS SUBROUTINE MEQSOL IS USE TO SOLVE THE MATRIX EQUATION AX=B OF 
N»IN THIS PARTICULAR CASFtN IS 4»ANO THE RESULT IS RETURNED TO 
THE MAIN PRGRAM BY THE LEBEU COMMON X(25).IF HIGHER ORDER MATRIX 
EQUATION is NEEDED TO BE SOLVE, THEN THE VALUE OF N HAVE TO BE 
CHANGE 


INTEGER M(25) 

DO 9 I=1,N 

M( I >a 1 

AMAX=A(I,1) 

DO 2 J«2,N 

IF<ABS(A( I, J ) )-ABS(AMAXn 2 . 2,1 
AMAX=A( I , J) 

M( I IsJ 
CONTINUE. 

IF(AMAX) 3,98,3 
00 4 J=1,N 
A(I*J)=AU , Jl/AMAX 
B( I ) = B< I )/AMAX 
DO 9 IP=1 , N 
IF( IP-ns,9,5 
MMM=M( I ) 

ZMULTaA(IP,MMM) 

DO 8 J=l,N 
IF{ J-MMM)7,6,7 
AflP, J)=0,0 
GO TO 8. „ . 

At IP,J)aA( IP, J)*“ZMULT4A{ I, J) 
CONTINUE 

et IP ) = B( IP >-ZMULT48( I ) 

CONTINUE 
DO 95 1=1 ,N 
XtM( I ) )=BC 1 ) 

RETURN 

WR1TE(6,99) .. 

FORMAT (12H NO SOLUTION) 

RETURN 

END 
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I 

ijj 


0001 

c 

c 

c 

c 

0 00'2 

0003 

0004 

0005 

0006 

0007 

0008 

0009_._ 

001 0 
001 1 
001 2 

0013 . • 

0014 
001 5 
0016 

0 017 

0018 

0019 

0020 

0021 . _ 


100 


SUBROUTINE MULTR ( N » X » Y » RMULT) 

THIS SUBROUTINE MULTR IS USE TO CALCULATE THF MULTIPLE 
CORRELATION COEFF R AND DENOTED THAT BY THE OUTPUT PARAMETER 
MULTR AND THE FOMULA IS R=SQRT( SUME OF P ( I t Y-) * B ( J 1 ) . 

WHERE B(J> IS THE STANDARDIZE REGRESSION COEFFANO THE 
UNSTANDARDIZE REGRESSION COEFF IS B=SDC J ) /S D( Y ) ( J ) 

REAL XdOOl.YUOO) 

EN=N 
SX=0. 0 
SY=0.0 
SXY=0.0 
SXSQsO.O 
SYSO=0 .0 
DO 100 1 = 1, N 
SX=SX+X( I ) 

SY=SY+Y( I ) 

SXY=SXY+X( I )4Y( n 
SXSO=SXSO+X( I )**2 
SYSQ=SYSQ+Y{ I )**2 
CONTINUE 

A=<EN*SXY-(SX1* (SY ) ) 
a=SQRT(EN»SXSQ- ( SX) **2) 

C=SORT(EN*'SYSO-(SY )*^2) 

RMULT=( A/(B*C) ) 

RETURN 

END 
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GRAPH 


0001 

0002 

0003 

0004 

0005 
0005 
0007 
0003 
0009 
COlO 
001 1 
0012 
0013 


SUBROUTINE GRAPH! Y) 
DIMENSION IP (.81 ) 

DATA lASTR/** • /♦IBLNK/* •/ 
J=ABS (Y ) 

IF( J-81 )500 .500»600 
500 DO 333 1=1.81 
IP(I)=IBLNK 
333 CONTINUE 

IP{ J > = IASTR 
.WRITE(6.S1) IP 
51 FORMAT! • +* ,38X. ' I • .01A 1 ) 
600 RETURN 
END 
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0009 


SUBROUTINE CORINE, R) 

COMMON/ApEAl/SX ,5V ,SXSQ, SVSQ, SXY 
C 

C IN THIS SUBROUTINE COR, THE INPUT ARGUMENTS ARE NE , SX ,SY ,SXSQ , SYSQ. 

C SXY; and the OUPUT argument is RJTHE COP, COEFF IS RETURN TO THE 

C MAIN PROGRAM BY R, 

C 

EN=FL0ATI NE) 

RNEW=EN4SXV-SX*SY 
DENEWX=EN*SXS0-SX>^+2 
0EN£WY=FNYSYSa-SY**2 
R=RNEW/SQRT (DENEWXYDENEWY) 

RETURN 

END 


LO 

o> 
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.C 

c 

c 


200 


SUBROUTINE STANDV ( N. Z . STDV ) 

THIS SUBROUTINE IS USE TO CALCULATE THE STANDARD DEVIATION 
OF A GIVEN SET OF DATA PTSAND DENOTED THAT BY THE OUTPUT 
PARAMETER STOV 
REAL Z(ICO) 

EN=N 
SZ=0 .0 
SZSQ=0.0 
DO 200 1=1. N 

SZ=SZ+Z( I ) 

SZSQ=SZSQ<-Zn )**2 
CONT INUE 

A=(EN*SZSQ-SZT*2 ) 

B-A/(EN+( EM-1 .0) ) 

STDV=SQRT(B) ’ 

RETURN 

END 
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F ^ £ R N A L S I G 

INTEGER VS ( 60 C ) . V V ( 60 0 ) , V ^ C 600 ) , V C ( 60 0 ) » VD^< AX ( 600 1 » VDM IN ( € 00 ) » 
ffiVP ( 600 ) t’STA TN , YEAR ,M0NTH.DAY,0^'AX^0M^tFReClP 
IN TEGER ^^e)♦^2(2>.^^{2)♦^4(3^.^5<3).N6{3),^7(3> 

INTEGER STN 1 , YRl ,RF 
REAL RVA )< ,R V IK 
DO 5CC IN=9.IC 

READ! IN,2G)M ♦^2t^2,^4 .^S ,^6,^7 
2C F0RMAT{eAl»2Al ,2A1 ,3A1 ,IX,3AJ «3Al .3X.3A1 ) 

I C N= 1 

EC REA0(IN*4C»END-499)£1A1NtYEAR»VCNTH,DAY'.C^AX»CVlN»PRECIPtRMAX» 

SIRM IN 

VS(N > = £TATN 

4C FORMA T(fe,3I2tlX,2I3,3X,I4,T14,2F3.0) 

VV ( N ) = VEAR 
VM I N ) = VaNTH 
VO ( N )=DA Y 
VP( N )=P«6CIP 

IFKDMAX.EQ.C) .ANO.I SIGNU.O.RMAXI.EC.-l.O)) CMAX=999 
VDMAM N )=OMA > 

IF( (DMIN.EQ.C) .AND.C £IG^(l•0,FM^)•EC.-l .0) )CV IN=999 
VDM INI-N )=DM I N 
M 1= W ( N ) 

M2= VD( N ) 

M 3= V S I K ) 

M4-=VS{ 1) 

IF { { (M 1 ,EQ. 12 ) .AND. ( ^2.EC.3n ) . CF. ( 0*3 .^E. V4 ) ) ) GO TO 200 
N = N+ 1 
GO TD EC 

2C0 CALL EORTNDI VS , VY , VM ,V0 .VO^AX .VOf' IN , N , ST N 1 . YR I , VP, P IB AR ) 

CALL ANP(STN1,YR1 ,VF ,BF,N) 

CLEAR ALL STCRAGF FCF NEXT STATICN 

DO 3CC 1=1. N 

V£( I )=C 

VY( I )= C 

VM( I )=C 

VD( I )=C 

vp( n=c 

3CC CCNTINLE 
GO TO 16 
4GG REIHINO IN 
5CC CONTINLE 
STOP 

CF6LG SLBCHK 
END 
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0033 
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0037 

0038 

0039 
OOAO 
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0042 
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C 

C 


THIS SLBROLTINF IS TC SCRTED THE 12 FCNThS CRSERVATIONS FOR 
EACH £^A^^C^4 

SL8R0L1INE SOR T MD ( V S » V Y i V K « VD i VO N A X « VDM h , N *S TM , Y R 1 , VP. P I BAR ) 
REAL 8M ,BM2 ,B V3 ,EV4 .BVS ,e N6 , E f'? » P V8 . BN9 » 8M 10 , BM 11 . BM I 2. 

4EN 1 .BN2 .8N3 ,P^4,E^£.B^e»B^7»e^a♦e^9«BM0»e^ll,P^12 
INTEGER OX 1( EC ) .DX2I EC) .0 X3 (SO ) .0X4 (50 ) .0X5(50 ) ,0X6 (SO ) , DX7<S0 ), 

♦ DXP(SC),DX9(EC),OXIC(5C),DX11(5C),DX12(SO),DM(50) ,DN2(50), 

4DN3( EC),ON4(5C),ONS(50) .0N6(50> ,DN7(50),CN8(50), ONQ (50 ), ON 10 { 50 ). 
4 ON 11(50). DN12(E0). 11. 12, 13, 14, 15, 16. 17, 18, 19.110.111, 112. YR I, YR2, 
4YR3,YR4.YPS,'NR6,>R7,VRP,VR9.YF10,YF11,YR12,STN1 ,STN2,STN3, STN4, 
4ETNE.£7NC,ST^7,ST^e,£T^9,S^M0,STMl,ST^l2 
integer \(S(eCC) ,VY(6C0) ,VM600) ,VC(600) ,VOVAX(600),VDMIN(600). 

4 JAN .FEB .MAR .APR , M A Y , ♦ J t L , AL'G , SE P ,CCT ,NCV «CEC, D1 (50 ), D2(50 ) , 

4O3(£C),D4(EC),O5(5C),D6<5C),D7(E0),Dfl(S0),C9(50>,Cl0(S0).DlH50>. 
SO 1 2( £0 ) .PI ( 50 ,P2( SC ) ,F3 (50) ,P4 (50 ) ,PS ( SO ) . P6 ( 50 ) . P7 ( SO > . P8 ( 50 ) . 
aPS(S0),P10(5C) .PIKSO) .912(50) ,VF(400 ) .FIEAP 
INTEGER Pi BAR ,P2eAR,93BAR,F4eAR,P5EAR.P6EAR.P7BAR, P8EAR, 
aP9BAn ,P 1 CHAR ,P I ISAR ,F1 2RAR 
I 1=C 
12=0 
13=0 
I4=C 
I5=C 
I6=C 
I7=C 

xe=c 

T9=C 
I 1C=0 


I 1 1 = C 
I 12=C 

DO ICO 1=1, N 
K=VM{ I > 

GO TOC 1.2, 3. 4, 5. £.7. 8. 9. 10. 11 .12) ,K 

1 11=1141 

P 1 ( I t )= VP{ I ) 

DX 1 ( II ) = VDMA >( I ) 

ON 1 ( I 1 )= VOMIN( I ) 

OK II )=VD( I) 
jan=vm( r ) 

VR 1=VV( I ) 

£TN1 = VS( I ) 

GO TO ICC 

2 12=1241 

P2( 12 ) = VP( I ) 

DX2( 12 ) = VDMA X I ) 

DN2( T2) = V0MIN( T > 

02( 12 ) = VD( I ) 

FE8=VM( I ) 

YR2=VV( I > 

STN2=V£( I ) 

GO TO ICC 

3 13=1341 

P3( 13 )=VP( I ) 

0X2( 13 )= VDMAX I > 

0N3( I3) = V0MIM 1) 

03(13 )=V0{ I > 
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00A4 

- 

MAR=VM( I ) 

0045 


>fR3= VY( I ) 

0046 


STN2=VS( 1 ) 

0047 


GO TO ICC 

0048 

4 

14=14+1 

0044 


P4( 14 )=VP{ I ) 

0050 


0>4( 14 > = VOMA X I ) 

OOS 1 


0N4( 14 ) = VDMIN( I ) 

00S2 


04 { 1 4 )= \/D( I ) 

0CS2 


APR= VM { I ) 

OOf 4 


VR 4= VV( I ) 

OOSS 


STN4= VS( I ) 

0056 


GO TO ICC 

0057 

c 

15=16+1 

0058 


P5( I 6 ) = VP( I > 

0059 


D>5( T5 >=V08A >( I ) 

0060 


ON6( IS ) = VDMIM I 1 

006 1 


D5( 16 ) = \D( I ) 

0062 


MA y=VM ( r ) 

0063 


VR 5= VT( I ) 

0064 


STN5=\iS( I ) 

0066 


GO TO ICC 

0066 

6 

16= I C+ 1 

W 0067 


P6( I 6 ) = VP( 1 ) 

1 006E 


DX6( 16 ) = VOyA >{ I > 

C 0069 


DN6( 16 > = VDM TN( I ) 

0070 


D6( I 6 ) = VD( I ) 

0 07 1 


JUN=VM( I ) 

007? 


YR6=Vy ( I > 

00 7 3 


STN6= VS{ I ) 

0074 


GO TO ICC 

00,7 5 

7 

17=17+1 

0076 


P7( I 7 ) = VP( I ) 

0077 


0>7{ 17) = VO VAX!) 

0078 


DN7 (17 > = V0^' TN( I ) 

0079 


07( 17 )=VD( I ) 

ooeo 


JUL=VM ( T ) 

008 1 


YR7=vy( r ) 

0062 


STN7=VS( I ) 

0C63 


GO TO ICC 

0084 

e 

18=18+ 1 

0085 


P6( 16 )=VP( I > 

C 066 


DX6( I e ) = VOVA X I ) 

0087 


DNe( 18 > = VDMIN( I ) 

0088 


AUG= I ) 

0089 


VR8=VM( I ) 

C090 


08( I € ) = VD( I ) 

009 1 


STN8= V£( I ) 

0092 


GO TG ICC 

C093 

9 

19=19+1 

0C94 


P9( I 9 >= VP( I ) 

0095 


D>9( 19 )= YDMA X I ) 

0096 


DN9 ( 19 ) = Vf)M IA( 1 ) 

0097 


09( 19 )=\D( I ) 

0098 


SEP= VM ( 1 > 

0099 


VR9= VX I ) 
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0 107 
0 loa 
0 109 
0 1 10 
0 111 
0 112 
0 113 
0 114 
Oils 
0 116 
0 117 
0 1 1 e 
0119 
0 120 
0 12 1 
“0 122 
M 0 123 
•^0 124 
0 125 
0 126 
0 127 
0 128 
0 129 
0 130 


0131 
0 132 
0 133 
0 134 
0 1 35 
0 1 36 
0 137 
0 138 
0 1 39 
0 140 
0 14 1 
0 14 2 
0 143 
0 14 4 
0 145 
0 146 
0 147 
0 146 
0 149 
0 150 
0 15 1 
0 152 
0153 


STN9=VS( I > 

GO TO ICC 

I C II C= I 1 C+ 1 

P1C(I1C) = VP( I) 

DX 1 C( I 1C ) = VOVA X I ) 

0N1C( I 1C )=VDVIN{ I ) 

0 1 C ( I I C ) = VD ( I) 

0C1=VM( I ) 

VR 1 C=VX I ) 

STN 1 C=VS( t ) 

GO TO ICC ' 

11 I I 1= I 1 1+ 1 

P 1 1( 1 1 1 ) = VP( I ) 

DX 1 1 ( 111 > = VDVA X T ) 

DN 1 l{ I 1 1 > = VONIM I ) 

D I l( I 1 1 ) = VD( I ) 

NOV=VM< I > 

YR 1 1 = VX I ) 

STNl l=VS( I ) 

GO TO ICC 

12 112=112+1 

P 1 2( I 12 ) = VP( I ) 

DX 12{ I 12 ) = VDVA X I ) 

DN 12( I 12 ) = VDNIM I ) 

012(1 12 ) = V0( I ) 

DEC= VM( I ) 

YR 12=VY( I ) 

STN 12=VS( I ) 

ICC CONTINUE 

WR ITE( e , EC ) 

50 FORMAT! 2CX »' THE VCNTHLY VAX, AND VIN TEVP TABLE*"/ 

32CX, ‘MONTH * *SX , *VA > TE VP • , 1 0 X , • V I h T E VP » » 1 OX , • P EEC IP • / 
a2CX,£2( •- * ) ) 

CALL OMEAN (OX I ,OM.Ol.JAN.YRl,STM,Il.EVl,eM) 

CALL OMEAN( DX2 tON2 ,D2 .FEB ,YR2 «STN2 t 12 ,e N2 , EN2 ) ’ 

CALL 0VEAN(0X3 ,DN3,D3,VAR,YR3,STN3.T3.BN3,EN3) 

CALL DMEAN (DX4 ,D^4,D4,AFfi,YR4,£T^4,14,e^4,E^4) 

CALL DVEAN(0X5.0N5.D5.VAY.YR5.STNS.T5.eN5.BN5) 

CALL DVEAN (0XC,0N6.D6.JLN.>H6.STN6,I6.EN6»EN6) 

CALL DMEAN< DX7,D^7,07,JLL^YR7,ST^7,I7,EV7»E^7) 

CALL DMEAN(DXe,DNe.De.ALG.Yfi8»£TN8.I8,ENe,6N8) 

CALL DMEAN (D >9 ,DN9 ,0 9 . SEP ,YR9 ,STN9,I9.eV9,EN9) 

CALL OVEAN( D X 1 C ,OMO ,01 0 .CCT ,YP10 ,STM0 , 1 10 .ev 10 .BN 10 > 
CALL 0MEAN(0Xll,DMI,0 11,NCV,YRll,STNll,Ill,e«ll,BNll) 
CALL DMEAN (0X12,0 M2 ,012,DEC,YP12,STM2,I12,evl2,BN12) 
CALL PBAR(P1 , I 1 ,PIBAP) 

CALL PEAR ( P2 , I 2 ,P2BAF ) 

CALL PBAR(P3 ,13 .R30AR) 

CALL PEAR(P4 , 14 ,P48AP) 

CALL PBAR(P5 , IS.FSBAF) 

CALL PEAR (P6 ,I 6 .F6HAPJ 
CALL PEAR (P7 , I 7 ,F7HAR ) 

CALL PE AR ( P8 ,I 8 ,F8RARJ 
CALL PE AR ( pc , I 9 ,F«3RAF) 

CALL PBAR (P I C , r 1 C ,F1 CRAR ) 

CALL PEAR< PI 1 ,I 1 1 ,P1 IRAR) 


O O- 

hsj) ftj) 




FCPTR/»N 

0 1 e4 
0 IS5 
0 1S6 
0 157 


0 158 
0 1 5S 
0 160 
0 t6 I 


IV G1 RELEASE 2.0 EQR7VD DATE = 77250 08/29/30 PAGE 

CALL PEARIP I 2 , I 12 ,F 1 2BAR ) 

V»RITE( 6.52>S1M *>R1 

52 FORMAT! 3)< S7A1! CN *'»I5»* IN 19', 12) 

VkRITEle.ei) JAN .HM 1 ,BM »F1 EAR , 
aFEB ,BM2 ,BN2 »P2BAR , 

3MAR ,BM3»BN3.P3B'AR , 

SAPR ,EM4 ,HNA .PAB AR . 
aMAT ,BM£ ,BNE .PSBAR , 
aJUN i8M€ »BN6 ,P6BAR , 
aJLL,SM7.RN7«P70AR. 

3AUG .8M5 ,HN6 tP8BAR « 
aSEP .BM9 .BN? ,P?BAn , 
aOCT.BM 1 C »BN 1 C »PI CBAR » 
aNOV.BM 11,BN11,P11BAR, 
aOEC ,BM 1 2 ,0N I 2 ,P12BAR 
51 FORMAT(2CX,I4,F12.2,10X.F10.2»10X,Ifl> 

RETURN 

DEBUG SLBCHK 
END 










A 



( END ') 
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FrRTR/»N 


000 I 
0002 

0003 

0004 

0005 

0006 
0007 
CC08 
0009 
00 10 
00 11 
00 12 
0013 
00 14 
00 1 E 


IV Gl RELEASF 2.C NAIN D/iTE = 77250 08/29/30 PAGE 

C THIS StBROUTTNE IS USE TC CALCULATE THE /i/CNTHLY WAX.MIN TEMP 

SLOROLT INE D^'EAN(D)^tO^*D,^<C^,TF^S7^♦^,^AXEAR, LINEAR) 

INTEGER DX{ 50 .DM SO) .MCN , 7R ,STh 
INTEGER 0(£O 
REAL MAXBAR ,KINHAR 
CALL TOTALIN iSMA> ,DX ,ICCLM) 

EN=N- ICCLNT 

IF( EN .eO.O .0 )EN= I . C 

MAXEAR^SRA >/EN 

CALL TC^AL(N,S^/I^,C^.KCL^T) 

EM=N-K0LNT 
IF(EM.EO.O.C)E^=1.0 
M INfiAR=SV I N/EV 
return 

DEOLG SL8CHK 
END 
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FCPTPAN IV Gl RELEASE 2.C 


7C7AL 


DATE = 77250 


08/29/30 


PAGE 


000 \ 
0002 
000 3 

0004 

0005 

0006 
0007 

ooca 

0C09 
0010 
00 1 1 
00 12 
0013 


StBROUTINE TDTALI N .StV ,DZ .KNT) 
INTEGER OZIEC) 

KN7=C 
SLM=C.C 
00 ICC K= I 

IF ( DZ( K ) .EQ .99? IGC 7C 200 
S1.M=SLM+0Z(K ) 

GO 70 ICC 
2CC KN7=KN7+1 
ICC CONTINUE 
RETURN 

DEBUG SU8CHK 
END 


w 

I 

to 

o 
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FCPTfiAN 


000 1 
0002 

0003 

0004 

0005 

0006 
0007 
OOOR 
0009 
00 10 
00 1 1 
00 t 2 
00 13 


00 14 
0015 
00 le 


IV Cl RELEASE 2,0 NAIN CATE = 77260 OF/29/30 PAGE 

C THIS SLBROLTINE IS LSE TC FIPO THE ANN FRECIF CF EACH STATION 

SUBROLTINE ANP ( S TAT N ♦ YEA P ,YP .EP . N) 

INTEGER STATN , YEAR ,VF{6Q0) .BP 
L=C 

IC0LN7=C 
DO ICC 1=1. N 

IF( VP( I ) .EG. SSSS) GC TC 200 
L=L+ VP( I ) 

GO TO ICC 

2CC ICOCNT= ICOLNT+ 1 
ICC CONTINUE 
BP=L 

t»RlTF(€.3C)N,ICCl.NT,STATN.YEAR,PF 
20 F0RMAT(3CX, MQTALA CF CB SE P V A T I C NS= • . 1 8 . / 
a2CX,»A0F MISSING CB SE R Y A T I C N = » . I 8 ,/ 

220X.‘FOR STATION « ,I e , • IN THE YEAP ‘.lA, 
a* THE ANNUAL PRECIP IS *,16) 

RETURN 

DEBUG SUBCHK 
END 
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FCFTR/JN IV Gl RELEASE 2,C 


PRAR 


SLBRDLTINE PRAR { PC .N ,PA V ) 

INTEGER PC( SC) »PAV 

I £=0 

OQ 1C I=1*N 

IFIPCC I ) ,EQ.SS£9) 6C TC 10 
I£=IS+PC< I > 

1C CDNTINLF 
EN=N 
PA V= I S 
RETURN 

DEBUG SUBCHK 
END 


w 

I 

N3 


000 1 
0002 
000 3 

0004 

0005 
0C06 
0007 
OOOE 
0009 
00 10 
00 1 1 
00 12 


DATE 


77250 


Oe/29/30 


PAGE 
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level 21.6 ( JUN 74 ) 


OS/360 FORTRAN H 


DATE 


n 

I 

00 


ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
f SN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 
ISN 


0002 

0003 

0004 
"0005 

0006 

0007 

OOOO 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
0017 
0016 

0019 

0020 


ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 


ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 

ISN 


0021 

0022 

0023 

0024 

0025 

0026 
0027 
0026 

0029 

0030 
'0031 

0032 

0033 

0034 


0035 

003b 

0037 

0038 

0039 

0040 

0041 

0042 


compiler options - NAMES MA IN » OPTsQO .L INECNT=55 » 51 ZE=0000K » 

SOURCE*EWCOIC»NOLIST,NOOECK»LOAn.MAP»NOEDIT*ID*XREF 

C JH£ AOV ANALYSIS 

C •»•»«••»«•*•» HY 

c ruddy h JEUN 

" C ' ■»»•»■«■«•«■«•» »«■ <M»’ «■«■«■*■»«•■»■» ■(}•*•»■» ■»«■»•»•»•«••»■ ■»■» W ■»<»• fl- 

REAL xnOO.lOO) »T(100) .GIIOO) 

INTEGER R(IOO) 

FTA«=6.93 
R(l)=12 
R ( 2)=12 
R (3) =12 
M=12 
N=3^ 

NG=3 
N0=0 

DO 100 1=1. NG 
NO=NO+R(I) ‘ 

CONTINUE 

DO 22 I=l ,M . 

READIS.iO) (X(I.J) .J=1.N) 

FORMAT (.3(F5,2.1Xn ' 

continue 

WRITE (fo, 98) 

FORMAT! ' 1 » .59X. »TA~LIEN(#54662) •/ 

043X. 'LOCKHEED' »5X,*U.S. AIR FORCE' ,5X.'*H0NG KOnG ROYAL'/ 
r«D43X.46{ •-' ) ) 

DO 11 1=1, M 

WRITE (6,9) (X(I ,J) ,J=1 ,N) 

FORMAT (43X'.F10.2".5X.'F12.2,5X,F 13.3) 

CONTINUE 
SUM=0,0 
SUM5Q=0 ,0 
DO 200 J=1,NG 
M=R { J) 

DO 200 1=1. M 
SUM=SUM.X (I , J) 

' ■■ SUMSO = Sl)MSQ + X tr. J)*<^2 
200 CONTINUE 

WRITE (6,10b) SUM.SUMSQ 

105 FORMAT (20 X, '5UM=» ,F10.2,2X, »5UMSQ=* .Fl0,2//) 

’C ■■ MO IS The tt OF OriSERVATIONS 
C MG IS THE ff OF GROUP 

ENO=NO 

SaSUM=SlJM««2 
' TSS= I SUMSQ-SOSUM/ENO) 

C TSS IS ThF TOTAL SUM SO OF OBSERVATIONS' 

DO 300 J=1,NG 
T (J)=0.0 
L=R ( J) 

no 300 1=1, L 
T(J)=T(J)+X(I,J) 


100 


10 

22 

98 


■"9 

11 


^8 

II 

I? 


a 


CO 



C-0 


ISN 0043 
ISN 0044 
■ISN 0045 
ISN 0046 
ISN 0047 
ISN 004ti 
ISN 0049 
ISN 0050 
ISN 0051 
ISN 0052 
ISN 0053 

I SN 0054 
ISN 0055 

ISN 0056 

ISN 0057 

ISN 005B 
ISN 0059 
ISN 0060 
ISN 0061 
ISN 0062 


ISN 0063 


300 CONTINUE 

DO 12 1=1, NG 
WRITE(6, 106)1, T(I) 

106 FORMAT(20X»»SUM OF GROUP ( • , 1 3 , • ) IS »»F10.2/) 

12 CONTINUE 
SUMGSO=0.0 
DO 400 1=1 ,NG 
e=R(I) 

SUMGS0=5UMGS0 + T ( I ) i»«2/E 
400 CONTINUE 

8SS= (SUMGSQ-SQSUM/ENO) 

C RSS IS THE SUMSO BETWEEN GROUP 
WSS=TSS-RSS 
NGM1=NG-1 

C NGMl IS THE DEGREE OF FREEDOM FOR THE WITHIN THE GROUP 
N0MN6=N0-NG 

C NOMNG IS THE DEGREE OF FREEDOM FOR WITHIN THE GROUP 
NOMi=NO-l 

C NOMl IS THE DEGREE OF FREEDOM FOT THE TOTAL OBSERVATIONS 
8MS=BSS/FL0AT (NGMl) 

WMS=WSS/FL OAT (NOMNG) 

FCAL=BMS/WMS 

WRITE(b,20) 

20 FORMAT (55X» ' ANALYSIS OF VARIANCE'/ 

+50X, 'DEGREE OF »,7X»'SUM OF »/ 

>-26X» 'SOURCE OF VARI AT ION '♦ SX ♦» FREEDOM », I OX* * SQUARE » »5X» 

♦ •MEAN SQUARE* ,5X* »F RATIO'/ - • 

♦26X,75('„' )/) 

WRITE (6, Hi) NGMl ♦HSS»BMS, 


ISN 0064 


ISN 0065 
ISN 0067 
ISN 0068__ 
ISN 0069 
ISN 0070 
ISN 0071 
ISN 0072 
ISN’ 007T” 


tMOMNG , WSS5 WMS iFC AL* 

♦ NOMi,TSS ■ 

21 FORMAT (25X, 'BETWEEN GROUP • , 1 2X , 1 3 , 9X ,F1 3 . 3 , 5X ,F 1 0 . 3/ 
+25X,' WITHIN GROUPS' ,12X,I3,9X,F13.3,5X,F10.3,5X,F10.3/ 

+?5Xt'T0TAL» »20X»I3,9X»F13.3/ . . 

♦2bX,77(»^»)/////) 

TF(FCAL.GF.FTAB) GO TO 101 
WRIT£(6,23) 

__23. FOBMAT(30X, 'SINCE. FCAL<FTA0, SO ACCEPI_THE HYPOTHESE!) 

GO TO 600 
IQl WRITE(6,24) 

24 FORMAT(30X, 'SINCE FCAL>FTAB,SO REJECT THE HYPOTHFSEM 
600 STOP 
END 



